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Abstract 

We compute the charge-density oscillations in a Tomonaga-Luttinger 
model with impurities for a short-range electron-electron potential. 
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In a recent paper [|IJ we introduced an extended version of the Non- 
local Thirring model (NLT) |J which includes a local interaction between 
the fermion current and a classical background field C^. Thus, the NLT 
interaction Lagrangian density 

Cint = \q 2 \ d 2 yJ^(x)V iM) (x,y)J^{y), (1) 
is now modified by adding 

c imp = v<pv - M( x )m (2) 

The name Ci mp refers to the fact that this density can be used to study 
the interaction between electrons and impurities ||. Indeed, if one uses ([[]) in 
order to describe the forward scattering of one- dimensional spinless electrons 
Bl 5 then the first (second) term in the rhs of (§) models the forward 
(backward) scattering between electrons and impurities. In particular, if we 
choose Ci(x) = together with 

C (x) = V5(x )8(x 1 -d)= M(x) (3) 

where V is a constant, our model coincides with the one recently considered in 
ref . to study Friedel charge-density oscillations in a Id Tomonaga-Luttinger 
liquid with an impurity located at x\ = d. 

We now start by considering the general case (C M (x) and M(x) arbitrary). 
Using a convenient representation of the functional delta and introducing the 
vector field (Please see ref. for details), the partition function of the 
model under consideration can be written as 

Z = J DA^det{i$ + t loPF + g£ + (p- M(x))e" 5[A] , (4) 

with 

S(A) = lj d 2 xd 2 yV^(x,y)A,(x)A,(y), (5) 
where Vr} is such that 

J d 2 z V^(z,x)V M (y,z) = 5 2 (x-y). (6) 
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As shown in jl| the partition function (^) is equivalent to a purely bosonic 
one, corresponding to a non-local Sine-Gordon model given by 



C-NLSG = ^{d^{x)) 2 + 

- f d 2 yd ll if{x)d ill) {x,y)d ll if(y) + C' imp , (7) 



with 



C! i mv = F^d^ifix) - a ° 2 cos (3cp(x) (8) 

F^{x) represents a couple of classical functions to be related to the C^'s 
and a (x) is to be related, of course, to M(x). Indeed, it is straightforward 
to show that the corresponding partition functions coincide provided that 
the following identities hold: 
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*F(PoV(i) + Pi v (o)) + P 2 MP 2 + d( )P 2 o + d w p 



(9) 



M(x)cos(2p F x 1 ) = ^p- (10) 



and 



>-i( n) r ( rj) P,C(p)-2e tluPu B(p) 



where 



2g 2 PvPyB(p)h(-p) = 
A(p) 

-/3 (11) 

P 2 + d( )Po + d (i)Pi 



a (p) = ^p 2 + \ [^(o 1 (p)p? + ^(i) 1 Ml > ( 12 ) 

^(P) = ^^(0) 1 »0 + ^) 1 »?], (13) 
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C(p) = [V (0) \p)-V a) \ P )}poPi, (14) 

A(p) = C 2 (p)-AA(p)B(p), (15) 

and 

D(p,x l ,y l )=J2(e tpx >-e™). (16) 

i 

As stated in the introductory paragraph, the path-integral identification 
depicted above allows to make contact with recent descriptions of Id elec- 
tronic systems in the presence of fixed (not randomly distributed) impurities 
|J. We shall now discuss this possibility. To be specific we shall consider the 
charge- density expectation value in a Tomonaga-Luttinger electronic liquid 
(V(i) = 0), with a non- magnetic impurity located al x\ = d, 

P( x ) = (^ dl( P ~ C0S ( 2 Pfx%) cos(2y/nip))c Bos (17) 

where 



Cbos = U d M x )) 2 + f d 2 ydMx) Ym v^ d My) 

-2i J d 2 yp F xx(V {0) (x, y)df)d 1 <f(y) 
+j-C (x)d 1 <p 

—Co(x) cos(2pi?xi) COS (2 -fwip). 

It is convenient to have an even Lagrangian density. This is easily achieved 
through a translation in the field ip. The result is 



£'bos = f (^(x)) 
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+ J d*yd lip (x)*^d lt p(y) (18) 
— Co(x) cos(2pFXi) cos(2 A /7ry9) 

where / is a classical function satisfying 

dof = (19) 

and 
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+ £ Jd 2 yV i0) (x,y)d!f 
+-^C + 2i PF Jd 2 y yi dfV {x,y) = 



(20) 



Thus we get 



p(x) 



—j=dif = - cos(2p F xi) cos(\/47r/(xi)) 
(cos(V47r^))r Sos , 



(21) 



where cos(2p^xi) cos(v / 47r/(^i)) is called the Friedel oscillation and A(x) = 



Since £'_b os is not exactly solvable, we shall use the well-known self consis- 
tent harmonic approximation (SCHA) which consists in replacing C'bos 



where m(x) is an arbitrary function, related to the impurity, to be varia- 
tionally determined. Also note that we have already considered the most 
interesting case of instantaneous potentials of the form V(p)(x,y) = 5(xo — 
y )U{x 1 - yx). 

Let us now compute the envelope A(x) in this approximation: 




by 



£ = ^*)) 2 + ^V + 




(22) 




(23) 



with 



J(x,y) 



2i^5 2 (x - y) 



(24) 
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We are interested in the case 



m(x) = m5(xi — d). (25) 

Going to momentum space and doing a translation in the field <p(p) — > 
4>{jp) + a(p), with a(p) a classical function, we find 

1 r d 2 p 



A(x) = exp -- J Pj(-p) a (p) = e*W> (26) 
where a(p) satisfies 



/ \ - J (P) [ d 2 q a(q) 
a{p) = 2R^p) ~ J j2^ miP ~ q) 2R(=¥) (27) 



with 



r (p) = -w+Pi 

I 7T 

m(k) = 2irm5(-k )e- lkld 

J{k) = -2i^e- ik - x (28) 

As we can see, the calculation of A(x) is reduced to finding the classical 
function a(p), that is, to solving the integral equation fl2"7p. This is an integral 
equation with degenerated kernel whose solution is fl3| 



with 



a(p) = 2^ + K ( Po )f(-p) (29) 



f(p) 
9{Qi) 



2tt 



1 - J-oo d<ng(qi)f{po,-qi) 
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It is straightforward to verify that a(p) can be written in a more compact 
way as 

W R(p) 1 7T + ml (p ,0) J V ; 

where 

I(Po,r)= r dgi COS{qir ] u( , . (32) 

Jo %8 + tf(i + *Qtel) 

with r = |d — 

So we finally obtained A(r) as a functional of the electron-electron po- 
tential: 

Mr) = exp -i r dp f/(p , 0) - . (33) 

7T J-oo \ 71 + mI(Po, 0) / 

This is the main formal result of our work. In order to check the validity 
of our formula we shall now consider a short-range interaction U{q{) = U = 
constant. It is convenient to express A(r) in the form 



A(r) = exp(T + W)(r) (34) 

where one has defined 

oo 



T — — — / d Po I( Po ,0) 



7T JO 



W{r) = ^L dp \ + ml(p ,0) (35) 



These integrals are easily performed yielding 



2 roc roo 

T = / dp J dq 1 



Wo ™J "pl + q!(l + ^) 
= Hj] X (36) 
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and 



7T J-oo 7T + ml(p , 0) 

1 mr _ 77&T . , , 

+ - e TT El {-—), (37) 
where E^x) is the exponential integral function, fi — > + and A — > oo are 

~2U 



infrared and ultraviolet cutoffs respectively, and X — J1 + 



At this point, in order to compare our results with previous computations 

A ^ X > A 2 



we consider the long distance regime defined by %f « 1, ^ >> 1. We get 



A(r) = C(A,/i,A)AV^ (38) 
which exactly coincides with [O] and 0] under the same regime. 
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